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Using the Bogoliubov-de Gennes equation, we study the quasi-particle spectra 
around a single vortex in a d-wave superconductor, where a magnetic field is 



. 

parallel to the c-axis. In the temperature region where the Ginzburg-Landau 

O 

theory is valid, we find that the local density of states preserves a circular 
symmetry when the symmetry of the superconducting order parameter is 
pure d-wave. It, however, exhibits a four-fold symmetry when the mixing 
of a s-wave component occurs. A peak with a large energy gap is found in 
the local density of states at the center of the vortex, which corresponds to 
the lowest bound state. Our results are consistent with a recent scanning 
tunneling microscopy experiment in an YBa2Cu307_,5 (YBCO) monocrystal. 
The breakdown of the Eilenberger theory in YBCO in particular and in the 
high-T c superconductors in general is discussed. 
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Typeset using REVT^X 



It is now well established that most of the properties of the high-T c superconductors are 
described in terms of a d-wave superconductor with possible exception of the electron- 
doped Nd2- 2; Ce : j:Cu04 0. Therefore it is of great interest to study vortex states in a 
d-wave superconductor j§-|[J. For example, in a magnetic field parallel to the c-axis, it 
was predicted that a square lattice of vortices tilted by ir/4 from the a-axis is the most 
stable for T < 0.8T C ||, which has been seen recently by a small angle neutron scattering 
and a scanning tunneling microscopy (STM) experiment || in YE^CusOr-^ (YBCO) 
monocrystals at low temperatures (T < 10K), although elongated in the b direction. One of 
the interpretations is that this distortion is due to the orthorhomibicity of YBCO, although 
there are alternative interpretations based on the (d+s) admixture PJlQl. 

Here we shall restrict ourselves to a single vortex in a d-wave superconductor, where a 
magnetic field is parallel to the c-axis. Early analyses of the quasi-particle spectra around 
a single vortex in a d-wave superconductor within the Eilenberger theory (a semi-classical 
theory of a superconductor) JTT] exhibits a surprising fourfold symmetry in the local density 
of states [|T^,|T^], which has not been seen experimentally These analyses rely on the 
Riccati equation which is obtained in ]r2| by simplifying the Eilenberger equation [TT[]. It 



is to be noted that, in the low-temperature region, where the Ginzburg-Landau (GL) 



theory is invalid, is also investigated by solving the Eilenberger equation self-consistently. 

We believe that the above Eilenberger theory, which is essentially a semi-classical ap- 
proximation, is invalid for the high-T c superconductors fl4 |. In order to clarify the point, 
we shall give an analysis of the experimental results. 

In YBCO II, no peaks are found in zero-bias tunneling conductance at the center of the 



vortex, in contrast to a s-wave superconductor NbSe2 |15| . Further, at the center of the vor- 
tex, there is a peak around E = 0.25A in the tunneling conductance, where A=280K is the 
superconducting order parameter of YBCO at T = OK. Then the most natural interpretation 
is that this corresponds to the lowest bound state for a vortex in a d-wave superconductor 
analogous to the one predicted by Caroli, de Gennes and Matricon Since the lowest 
bound-state energy is given approximately by Eq = A/n^pp^^ 1 [12,|19|, where pf and £ are 



the Fermi momentum and the coherence length respectively, this implies pf£ ~ 1 in YBCO. 
This is also consistent with the chemical potential of YBCO deduced from the analysis of 
the spin gap seen in an inelastic neutron scattering experiment from monocrystals of YBCO 



Therefore in order to describe the quasi-particle spectra around a single vortex in high-T c 
superconductors, it is crucial to study the Bogoliubov-de Gennes equation which is given by 

{~Y~ ( V - zeA(x)) 2 - /i}m„(x; k) + A(x; fc)u n (x; k) 
= e„(/c)M n (x; k), 

-{-W~( V + ^A(x)) 2 - AiWx; k) + A*(x; fc)« n (x; k) 

= e n (k)v n {x;k), (1) 

where «„(x; k) and f n (x; k) are the quasiparticle amplitudes, A(x; k) is the pair potential, 
A(x) is the vector potential which is neglected assuming H<^H c2 ||16|| , and fi is the chemical 
potential which is identified with the Fermi energy. In this paper, we consider two types of 
the pair potential A(x; k), which are 'pure d-wave' 



A(x; k) = A tanh(r/£)cos(20(A;))e 10 , (2) 

where r, <p and 9{k) are defined by x = (r cos(0), r sin(0)), (k x ,k y ) = (cos(9(k)),sin(9(k))), 
A is a real constant and £ is the coherence length of the superconductor; and '(d+s)-wave' 
2| 



A(x; k) = Atanh(r/0 cos(20(fc))e^ + b {r/C)e~ iq 

for < r<£, 
A(x; k) = (Atanh(r/0 + a (£/r) 2 ) cos(20(fc))e^ 

+(£Af- 

cos(20(A;))(a_ 1 e- 3 ^ + a 1 e 5 ^) 
(6oe-^ + 6ie 3 ^) 



for £ < r. (3) 
The pair potentials are obtained by the GL theory for an anisotropic superconductor 



T2| , |20H and applicable at not too low temperatures (an estimate of the temperature region 
where the GL theory is valid is [0.5T C , T c ] for the s-wave superconductors, when one consider 
the quasi-particle spectra around a single vortex |19| ) . 



In writing ([!]), we neglected the non-commutability between k and x [21| |22 |. The correc 



tion is 0{1/pf^) and irrelevant at least in the study of systems with a long coherence length 
e.g. a superconducting phase of a heavy-fermion system (p F £ ~ 10) and 3 He superfluidity 
(PfC, ~ 100), but may have a serious influence in the study of the high-T c superconductors, 
where pp£ ~ 1 as is discussed above. But our results seem to be consistent with experi- 
mental results of YBCO as will be discussed below. Here it is to be noted that, in contrast 
to a usual semi-classical theory of a superconductor (the Eilenberger theory), our method 
takes into account quantization of the angular momentum. In fact, as we shall see, the 
quasi-particle spectra is totally different from ones obtained semi-classically [|l2l , |l3" . 



In order to solve the Bogoliubov-de Gennes equation (|l|) numerically, it is convenient to 
expand the quasi-particle amplitudes w n (x; k) and t> n (x; k) as 

oo oo 

l=— ooj=l 
oo oo 

v n (x;k) = J2 v n,i,j(k) ^|/_i|(r)exp(i(Z- 1)0). (4) 

/=— ooj=l 

Here it>i V (x) = 75=^-r- , — rJ u (ai U x/R) {J v (x) is the Bessel function), «, „ is the j th 
positive zero point of J v {x) and R is the radius of the system. In our numerical calculation, 
the number of the basis is about 80000 for pure d-wave and about 4800 for (d+s)-wave. It is 
sufficient to discuss physical properties around a single vortex in a d-wave superconductor. 

The quantity of interest for comparison to STM experiments is a local density of states 
in a superconductor 

AT(£,x)=p KM)| 2 5(E-e n (k)) 

k,n 

+ k(x,A;)| 2 8{E + e n (k))}. (5) 



The local densities of states N(E, x) for pure d-wave are shown in Figs. 1-4 for r/£ =0.0, 
1.0, 2.0 and 10.0, respectively. The parameters are chosen as ppC, = 1.41 and 2m£ 2 A = 2.82, 
from the experimental condition ||. We also set = 40. One can see peaks in the local 
densities of states corresponding to bound states around a single vortex, which are localized 
near the vortex center. 

For r/£ = 0.0 and 1.0 (Figs. 1 and 2), a peak with a large energy gap, which corresponds 
to the lowest bound state, appears in the local density of states. The peak has a width due 
to the internal degree of freedom in the k space. The result is consistent with a recent STM 
experiment ||. For r/£ = 10.0 (Fig. 4), N(E,x) : s around a vortex and without a vortex 
are almost identical. This indicates that the bound state is localized within this distance. 
We believe that the large energy gap is closely related to the quantization of the angular 
momentum. We note that there is a clear asymmetry in the local density of states at E and 



-E due to the formation of the bound states 23 



The integrated local densities of states near the zero energy are shown 

in Figs. 5 and 6 for pure d-wave and (d+s)-wave respectively. The parameters are chosen 
as p F £ = 1.33, R/i = 30, 2m£ 2 A = 2.82, 2m£ 2 a m = 0.05 (m = -1,0, 1) and 2m£ 2 b m = 
0.05 (m = 0, 1). The choice of high energy cutoff does not change the results quantitatively. 
Note that the local density of states preserves a circular symmetry for pure d-wave. On the 
other hand, it exhibits a clear four-fold symmetry when the mixing of a s-wave component 
occurs ( see Fig. 6 ). 

In conclusion, we have investigated the quasi-particle spectra around a single vortex in 
a d-wave superconductor by solving the Bogoliubov-de Gennes equation numerically. The 
local densities of states are consistent with a recent STM experiment || . The result shows a 
peak with a large energy gap in the local density of states at the center of the vortex, which 
corresponds to the lowest bound state. We find that the local density of states preserves 
a circular symmetry when the symmetry of the superconducting order parameter is pure 
d-wave. It, however, exhibits a four-fold symmetry when the mixing of a s-wave component 
occurs. However it is possible that a four-fold symmetry of a local density of states without 



a mixing of a s-wave component |l^ , [T3|| appears due to the noncommutability between k and 
x. A fully quantum-mechanical treatment of the Bogoliubov-de Gennes equation is left as 
a future problem. Also a self-consistent treatment is needed in low temperatures, where 
the GL theory is invalid. Moreover, a treatment of a vortex lattice is essential when one 
investigates the thermodynamic quantities. 

One of the authors (K.M.) thank 0. Fischer and Ch. Renner for illuminating discussions 
on their STM results. This work is in part supported by National Science Foundation under 
grant number DMR92-18371 and DMR95-31720 and by Grant-in- Aid for Scientific Research 
from the Ministry of Education, Science and Culture, Japan. 
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FIGURES 



Fig. 1 The local density of states N(E, x) around a vortex (o) and that without a vortex 
(+), where r/£ = 0.0. The rescaled energy is defined by E/A. 

Fig. 2 The local density of states N(E, x) around a vortex (o) and that without a vortex 
(+), where r/£ = 1.0. The rescaled energy is defined by E/A. 

Fig. 3 The local density of states N(E, x) around a vortex (o) and that without a vortex 
(+), where r/£ = 2.0. The rescaled energy is defined by E/A. 

Fig. 4 The local density of states N(E, x) around a vortex (o) and that without a vortex 
(+), where r/£ = 10.0. The rescaled energy is defined by E/A. Weak anisotropy appears 
due to boundary effect. 

Fig. 5 The integrated local density of states J °' 17A dE N(E, x) for pure d-wave (left) and 
that for (d+s)-wave (right), where = 3ir/8. The lines are guide for eyes. 

Fig. 6 The integrated local density of states Jq' 17A dE N(E,x) for pure d-wave (o) and 
that for (d+s)-wave (+), where r/£ = 3.0. The lines are guide for eyes. 
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